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We explore the phase structure of two-color and two- flavor QCD in the space of the quark chemical 
potential /iq and the isospin chemical potential yci. Using a mean-field model we calculate the chiral 
and diquark condensates, a and A, self-consistently. In weak coupling and in the chiral limit, we 
confirm the interval of the isospin chemical potential, 0.7IAo < fix < 0.75Ao, in which a single 
plane-wave Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) phase is favored over isotropic superfluidity 
and normal quark matter. The LOFF window becomes slightly wider at high density. For stronger 
coupling with nonzero quark mass, which is relevant to currently available numerical simulations in 
lattice two-color QCD, the single plane-wave LOFF phase appears only at sufficiently high density. 
The prediction obtained for the LOFF region could be tested with lattice since we can prove that the 
present system is free from the fermion sign problem. We draw the energy landscape on which local 
minima corresponding to the isotropic superfluid phase and the LOFF phase and a local maximum 
corresponding to the gapless phase are manifest. Our results clearly illustrate the path from the the 
unstable gapless phase down to the LOFF phase. 

PACS numbers: 12.38.Mh, 26.60.+C 



I. INTRODUCTION 

It is a longstanding problem to uncover the phase 
structure of dense nuclear and quark matter in the low 
temperature and high baryon density region because of 
its complexity. From the academic point of view, our 
curiosity urges us to imagine what an extreme state of 
cold quark matter at asymptotic high density is like. In 
fact, at density far larger than the strange quark mass 
M s but still smaller than the charm, bottom, and top 
quark masses, there has been established a consensus 
that quark matter takes on color superconductivity in 
a color-flavor locked (CFL) manner [l], EJ. Then, what 
comes next as the density goes down? This is an im- 
portant question because, if quark matter appeared in 
neutron star cores, its state would be strongly affected 
by M s . One plausible candidate was considered to be 
a gapless color superconducting phase [H, which is a 
QCD version of the Sarma phase [1, Q partially stabi- 
lized by neutrality constraints. Specifically the gapless 
CFL (gCFL) phase [6| was expected in quark matter in 
the intermediate density region. It turns out, however, 
that the gapless phase is unlikely to exist in such quark 
matter because of the chromomagnetic instability that 
develops at sufficiently low temperatures 0, H H EH • At 
higher temperatures, a it-quark superconducting (uSC) 
phase is predicted to occur as a remnant of the gapless 
phase [H EH, and the existence of the doubly critical 
point facing both the uSC phase and a d-quark super- 
conducting (dSC) phase [T3| seems robust [la. fl4|. 

Interestingly, the chromomagnetic instability in the 
gapless phase tends toward spontaneous generation of 
a total momentum 2q carried by each Cooper pair [H 
El EBj- [In addition to the chromomagnetic instabil- 
ity, an instability appears with respect to inhomogeneous 
fluctuations in the gap magnitude, leading the gapless 
state to a BCS-normal phase separation or a BCS-normal 



mixed phase [17|. See also Ref. [l8|.] In general, an 
inhomogeneous superconducting phase characterized by 
pairing with nonvanishing 2q is referred to as the LOFF 
phase named after Larkin-Ovchinnikov [l9| and Fulde- 
Ferrell 20]. In the present context, it is possible to de- 
scribe a state resulting from the same instability in sev- 
eral different ways: a state with current generation of 
collective excitations [2l], [22| , a state with gluon conden- 
sation 23], and a colored LOFF state H 16| . They are 
all equivalent algebraically. In terms of the colored LOFF 
state, the chromomagnetic instability is to be regarded 
as an instability with respect to the spatially oscillating 
phase factor, e lX q r with the Gell-Mann matrices A Q , 
of the pairing gap matrix in color space. In the present 
paper which focuses on a two-color theory, we shall gener- 
alize the notion of chromomagnetic instability to a phase 
instability that makes sense not only in a charged super- 
conductor but also in a neutral superfluid. 

In the LOFF phase the translational and rotational 
symmetries are spontaneously broken by q. The sin- 
gle plane-wave LOFF phase is only the simplest etude 
and in general a complicated crystal structure should 
emerge. In the context of QCD physics [H El El El, 
the plane-wave LOFF phase, crystal structure, stability, 
and its physical property have been examined mainly 
by means of the high-density effective theory (HDET) 
and the Ginzburg-Landau (GL) expansion in terms of 
the pairing gap [H El El Elj . For the moment it is 
an urgent problem to clarify the energetically favorable 
structure of the three-flavor LOFF phase. Whereas the 
phase instability guarantees existence of the LOFF phase 
with a lower energy than the gCFL phase even within the 
single plane-wave ansatz, the HDET and GL approxima- 
tions do not allow us to identify the most favorable LOFF 
phase in full details. 

We shall here revisit the LOFF phase in two-color and 
two-flavor QCD in the presence of both the quark chem- 
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ical potential /i q and the i sosp in chemical potential n\. 
It was pointed out in Ref. [32] that the interval of /ii in 
which the LOFF phase occurs exists in such a system at 
least in weak coupling, while it is nontrivial whether the 
LOFF window should survive or not at stronger coupling 
with a finite quark mass introduced. We will make use 
of a mean-field model to address this issue by following 
a line of argument of Ref. [33j, which looks successful 
in reproducing the numerical results of lattice two-color 

qcd urn, is us Him. 

Although it is beyond our current scope, we expect 
that the lattice two-color QCD could observe a signature 
of the LOFF phase numerically. We adopt two sets of the 
model parameters, one of which is relevant to currently 
available lattice simulations in which the diquark con- 
densate has been measured at /i q ^ but fxi = 0. This 
work is a first step toward identifying the LOFF state on 
lattice. 

We can also mention that it is instructive to shed light 
on two-color QCD as a mimic of real QCD in which the 
LOFF state is a natural consequence of the phase insta- 
bility in the gapless state. In the context of real QCD, 
we know the direction in which the unstable gapless state 
goes and consider the LOFF state as a natural replace- 
ment of the unstable state, but it remains to be surveyed 
how one gives way to the other. At this point, it is advan- 
tageous to switch to two-color QCD. In fact, simplicity 
inherent in two-color QCD enables us to picture the en- 
ergy landscape including various states. By doing so, we 
can get a feeling that we are heading for the right way 
dictated by real QCD. 

This paper is organized as follows. In Sec. [Til we men- 
tion unique features of two-color QCD, our mean-field 
model, approximations to be made, and the twofold pa- 
rameter choice at weak and intermediate coupling. Sec- 
tion lllll is composed of four subsections presenting numer- 
ical results from the mean-field model for the isotropic 
superfluid phase, the unstable gapless phase, the LOFF 
phase, and the energy landscape, respectively. Our main 
results are summarized in Fig. O for weak coupling with 
massless quarks and in Fig. [5] for intermediate coupling 
accessible in the lattice setting. We also plot the free en- 
ergy in the space of the diquark condensate and the pair 
momentum in Figs. and [5] for the intermediate coupling 
case. Section HVl is devoted to our conclusions and future 
perspectives. 



II. SETUP 

In this section, we will briefly summarize a chiral sym- 
metry breaking pattern inherent in two-color QCD, dis- 
cuss the absence of the fermion sign problem, and then 
describe a model for two-color matter with nonzero /i q 
and j«i along with some approximations in calculating the 
thermodynamic potential. We will emphasize the benefit 
from two-color nature that the model prediction is more 
robust than in three-color case. 



A. Symmetry Breaking Pattern 

Two-color QCD has a peculiar feature of chiral symme- 
try. When quarks are massless and the chemical potential 
is zero, chiral symmetry in two-color and two-flavor QCD 
is augmented from the standard one SUl(2) x SUr(2) to 
SU(4). This is because a doublet and an anti-doublet of 
the SU(2) group are indistinguishable, which makes left- 
handed quarks and right-handed antiquarks belong to 
the same group multiplet and hence doubles the basis. In 
fact this extra symmetry, which is often referred to as the 
Pauli-Giirsey symmetry [40j], amounts to rotational sym- 
metry among the chiral and diquark condensates. The 
spontaneous chiral symmetry breaking pattern is thus 
SU(4) — + Sp(4) in the presence of diquark condensation. 
Once a nonzero quark chemical potential sets in, i.e., 
/i q 7^ 0, the extended symmetry is reduced to standard 
global symmetry as SU(4) -> SU L (2) xSU R (2) xU B (l). If 
quarks are massive, the remaining symmetry is as usual 
SUy(2) x Ub(1), which spontaneously breaks down to 
Sp(2) once the quark chemical potential exceeds the mass 
of the lowest-lying excitations carrying nontrivial baryon 
number. Because the Ub(1) symmetry is spontaneously 
broken by the diquark condensation, the system is a 
superfluid accompanied by massless Nambu-Goldstonc 
bosons. The gauge symmetry is intact since the diquark 
condensate is gauge invariant in this case. 

It is obvious that no chromomagnetic instability can 
happen in any kind of superfluid state of this system be- 
cause no finite Meissner screening mass arises without 
color symmetry breaking. In other words, if we interpret 
it as the phase instability as we mentioned above, the 
color singlet diquark condensate cannot have any spa- 
tially oscillating phase in color space and thus cannot 
lead to any instability involving color degrees of freedom. 
In terms of the phase instability, on the other hand, we 
can understand that fluctuations in the Ub(1) phase in 
the gapless superfluid state could result in instability. It 
should be mentioned that the modes in the Ub(1) phase 
itself in a superfluid correspond to the massless Nambu- 
Goldstone bosons, and the instability would not emerge 
in the phase factor e l2q r but in the wave-vector q of the 
phase, i.e., the current of the Nambu-Goldstone boson in 
a physics language (see Refs. [2l|, [22| for details). 



B. Fermion Sign Problem 

In this subsection, we briefly summarize the sign prob- 
lem of the fermion determinant [4lj and make sure that 
two-color QCD escapes from it for any number of quark 
flavors. The Dirac operator in question reads 

M(flq) =7/^+74/^ + ™ (1) 

per one quark flavor in Euclidean space, where de- 
notes the covariant derivative containing gauge fields. In 
the presence of the isospin chemical potential fii with u 
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and d flavors, the Dirac operator is then a direct sum of 
two flavor sectors: A4(u u ) © For the moment 

we shall focus on the one-flavor A4(/i q ) to discuss the 
sign problem, since det A4(p q ) > is sufficient to claim 
det[M{u u ) © M(fid)} = detM(u u ) ■ detM(u d ) > 0. 

Because the Euclidean gamma matrices 7^'s are her- 
mitian by convention, the eigenvalue A„ of anti-hermitian 
1pP», i.e., 

= X n ^ n , (2) 

is pure imaginary. The eigenstate ^tpn has an eigenvalue 
— A„, which in turn equals to A* . The Dirac determinant 
is, therefore, a product of all the paired eigenvalues A„ + 
to and A* +to, which is non-negative because |A„ + to| 2 > 
0. 

When /i q is nonzero, 74/iq is hermitian and thus the 
eigenvalue A„ defined by 

(7 M L> M + 74/iq)V'r 1 = AnV'n (3) 

is no longer pure imaginary but complex. Here again, 
75-0n has an eigenvalue — A„, but it is different from A* 
in this case. Therefore, the Dirac determinant, given by 
n„(A„ + m)(— A„ + to), is not necessarily non-negative, 
and when it is negative for some gauge configurations, 
the sign problem occurs. 

Two-color QCD is unique in the sense that one can 
find another eigenstate a^C -1 757/)* with an eigenvalue 
A* , where 02 is the second Pauli matrix in color space 
and C represents the charge conjugation. Obviously, 
the eigenstate multiplied by 75 has an eigenvalue —A*. 
Consequently, the eigenvalues of the Dirac operator al- 
ways constitute a quartet: A„ + m, — A„ + m, A* + to, 
and — A* +to, leading to non-negative Dirac determinant 
through the relation |A„ + to| 2 | — A„ + m\ 2 > 0. We 
note that our argument holds even with odd number of 
quark flavors, although in this case the sign problem was 
supposed to occur [321] . 

We can thus conclude that two-color QCD has no sign 
problem at finite density regardless of the number of quark 
flavors. This implies that the Monte-Carlo simulation for 
two-color QCD with two flavors is feasible in the presence 
of both /i q and /xi, which would help us elucidate the 
phase structure. 

C. Model for Two-Color Quark Matter 

We turn to a model for two-color and two-flavor quark 
matter within the BCS ansatz that two particles with op- 
posite momenta p and —p pair with a short-range inter- 
action. We may then assume that only two mean-fields, 
namely, the chiral and diquark condensates, 

<j = G(^), A = Ge^e ab (^ a zC l5 ^) , (4) 

are predominant in the region < [i q . We denote the 
color and flavor indices by a, b and i, j, respectively. 



Here, T represents the transposition in the Dirac index, 
and C is the charge conjugate matrix to make the di- 
quark condensate Lorentz invariant. The dimensional co- 
efficient G, which controls the strength of the four-fermi 
coupling constant, is common to the chiral and diquark 
condensates by virtue of the Pauli-Giirsey symmetry. In 
a rough but intuitively more understandable notation, 
o~ ~ (uu) + (dd) and A ~ (ud). Note that when /zi > /i q , 
the predominant diquark condensate (ud) gives way to 
(ud). In this work, however, we limit our discussion solely 
to the region fj,j < fj, q . 

In the LOFF phase where a u-quark with the momen- 
tum q + p and a d-quark with the momentum q—p form 
a Cooper pair with the total momentum 2q, the diquark 
condensate has a spatially oscillating phase as 

A -> Ae l2qr . (5) 

It should be noted that we neglect any spin-one conden- 
sate such as (uu), (dd), and (dd) entirely in this work. 
It generally coexists with the spin-zero condensate (ud) 
in the LOFF state, but is known to be smaller by one 
order of magnitude than the spin-zero one [24[ ■ We re- 
mark that the spin-one condensates could be relevant for 
even larger Fermi surface separation between u and d 
quarks [32j. 

Then, the thermodynamic potential with the mean- 
field condensates a and A and the pair momentum q can 
be expressed as 

0(cr, A,q;m ,T, jU q ,/ii) 

2 J (2^) 3 ^\ 2 L J J 



Here A is the cut-off parameter and too is the current 
quark mass. After the p-integration, the thermodynamic 
potential no longer depends on the direction of the three- 
vector q but becomes a function of its magnitude q = \q\ 
alone. 

Expression ([6]) contains the sum of the quasi-quark en- 
ergies ej(p), which correspond to the 32 eigenvalues of the 
32 x 32 quark Hamiltonian matrix with two colors, two 
flavors, two spins, quark-antiquark, and Nambu-Gor'kov 
doubling. The remaining part is the energy shift asso- 
ciated with the mean-field approximation. The conden- 
sates, a and A, share the same coupling G as implied by 
the Pauli-Giirsey symmetry; in the chiral limit too = at 
A'q = Mi — 0, the thermodynamic potential f2 should be 
reduced to be a function of a 2 + A 2 . This feature exem- 
plifies a great advantage of two-color QCD over various 
model studies of real QCD in which quantities affected 
by chiral dynamics strongly depend on the model param- 
eters. 

When A = ui = 0, all quarks have the same con- 
stituent mass, M = too — cr, leading to the quark and 
antiquark energies, =F Mq with = \/p 2 + M 2 . 
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The effect of nonzero /Lti is to shift fi q by and — fii 
for u and d quarks, that is, 

M« = Mq + Ml I Md = Mq ~ Ml 1 (7) 

which results in crossing of the energy levels. In the pres- 
ence of the diquark condensate A, the crossing energy 
levels mix together and a level repulsion occurs. We will 
simply approximate the energy levels by mixing between 
the w-quark energy £(q + p) — /x u and the d-quark energy 
£(9 — P) ~ Md- Of course we can also swap q + p and 
q — p but the integration over p washes out the differ- 
ence eventually. The eigenvalues of the following 2x2 
matrix, 



£(<?+p)-M« A 

A -£(q -p) + /i d 



give the quasi-quark energy dispersion relations, 



e p (P. <?) = 



6£(P, q) ± Mi + V [£(p, 9)-Mj + A 2 



(8) 



(9) 



with 



^(p,<z)=£(<Z+p)-£(<?-p), (10) 

£(p,9) = M£(9 + p) + £(9-p)] ■ (11) 

Note that there are eightfold degeneracies for e+ (p) and 
e~(p), respectively, which implies sixteen energy levels 
in total. We should remark that these energy dispersion 
relations are slightly different from those used in Ref. [24| 
because the gap energy A is constant in our approxima- 
tion instead of depending on the relative angle of q + p 
and q — p. It is straightforward to introduce any q de- 
pendence in A. We will, however, stick to a constant A, 
partly because we will confirm later that a constant A 
is enough to reproduce results quantitatively consistent 
with Ref. and partly because no reliable ansatz is in 
hand a priori except when M is zero. 

As for the antiquark contribution, dropping A is a 
good approximation as long as A <C M + /i q . Such 
an approximation, however, would apparently ruin the 
Pauli-Giirsey symmetry which should be present in the 
limit of mo = /i q = /ii = q = 0. Hence, we shall keep A 
as well as a also for the antiquark energy levels by adopt- 
ing the same form of the energy dispersion relations as 
the quark contribution with — /i q replaced by +/tt q : 



: (P, <?) 



y^(p,q)±m + J[^(p 1 q)+H q \ +A^ 



(12) 



with eightfold degeneracies again, although difference of 

et(p,0) from l£(<Z + P)+M«l and |£(g-p)+Md| is merely 
negligible. We do not consider the mixing between quarks 
and antiquarks which would take place with nonvanishing 
(ud) when /ii > /i q . 



Using the thermodynamic potential as specified above, 
we will solve the following equations: 



on 



= 



dA 



= 0, 



dq 



= 0, 



(13) 



to obtain the self-consistent values of A, a, and q. Now 
that we have come by all the necessary formulae, we will 
proceed to numerics next. 



III. RESULTS 

In this section we will present the numerical solutions 
to Eq. (fl"3|) obtained for the isotropic superfluid phase, 
the gapless phase, and the LOFF phase. The gapless 
phase is unstable, but still we will investigate the gapless 
solution to Eq. (fT3|) for the later purpose of surveying the 
energy landscape. 

For numerical evaluations we first have to fix the model 
parameters, namely, the current quark mass mo, the cut- 
off A, and the four-fermi coupling constant G. We do 
not need to specify A because we will present all the 
dimensional quantities in unit of A. Instead of G we 
will use the diquark condensate at fj,i = denoted by 
Ao(/Lt q ) to specify the interaction strength; a larger Ao 
means a stronger coupling. In this paper we will take two 
parameter choices: 



Parameter I: m = , A (^ q = 0.5A) = 0.05A, (14) 



and 



Parameter II: m = 0.025A , A (/i q = 0.5A) = 0.5A . 

(15) 

Hereafter we will refer to Parameter I as the weak cou- 
pling parameter set since it will allow us to confirm that 
our approach recovers the known properties of the LOFF 
phase at weak coupling. For Parameter II, which will be 
denoted by the intermediate coupling parameter set, the 
choice of mo is motivated by available lattice simulations 
of two-color QCD. As we can see in Fig. [H for Parameter 
II, the chiral and diquark condensates amount to a com- 
parable magnitude. In the absence of the sign problem 
of the Dirac determinant even with fj,\ introduced as dis- 
cussed in Sec. Ill B| our results at intermediate coupling 
could be readily compared with what we would observe 
from lattice simulations. 



A. Superfluid Phase 

We show in Fig. [T] the behavior of the chiral and di- 
quark condensates as a function of the quark chemical 
potential fj, q at /ii = q = 0. The condensates at weak 
coupling with Parameter I (|14p are exponentially sup- 
pressed for small /i q . This suppression is a sharp contrast 
to the results in the strong coupling limit [42l . [43j | . In par- 
ticular a keeps vanishing entirely because of mo = as 
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Quark Chemical Potential ju q [A] Isospin Chemical Potential n x [A] 



FIG. 1: Chiral and diquark condensates as a function of iz q 
at Hi = q = 0, given in unit of A. The quantities with the 
subscript "I" and "II" are the results obtained at weak and 
intermediate coupling, respectively. 



FIG. 2: Superfluid and gapless solutions at fiq — 0.5A. The 
first-order phase transition from superfluid to normal quark 
matter is indicated by the dotted line. 



first noted in Ref. [42| ■ From Fig. Q] we can make sure 
that Ai(^ q = 0.5A) = 0.05A is satisfied. In the case of in- 
termediate coupling with Parameter II (fT5|) , on the other 
hand, the diquark condensate appears above the thresh- 
old fiq = 0.185 A. We can confirm that An(^ q = 0.5A) = 
0.5A as it should. The non-monotonic behavior of A at 
f-q ^ 0.7 is caused by the saturation effect: The den- 
sity cannot grow unboundedly because the phase space 
is limited by the presence of A, which is observed in lat- 
tice simulations as well as in model studies [H, |t| . We 
remark that our An and |cth| are to be compared with 
the results given in Figs. 1 and 2 of Ref. [33]. 



B. Gapless Phase 

The isospin chemical potential fii plays a role in ex- 
erting stress onto quarks to tear the Cooper pair apart. 
As long as jUi < A, the isotropic superfluid phase is rigid 
against such stress. This can be easily understood from 
Eq. §§§ with q — substituted; if /ii < A and q = 0, then 
e+ + £p = (£ — fj, q ) 2 + A 2 , which has no dependence 
on Hi, and thus A and a stay constant regardless of fii- 

Once /ii becomes larger than A, e~ decreases down to 
zero at p ~ /i q , that is, the energy dispersion relation be- 
come gapless. This gapless superfluid solution is similar 
to the gapless 2SC (g2SC) phase known from QCD. In 
contrast to the g2SC case, we do not impose neutrality 
on the present system. In principle, we can consider elec- 
tric neutrality by turning the electromagnetic interaction 
on and putting electrons into the system. We will not do 
so, however, partly because the resultant stability of the 
gapless state with respect to homogeneous change in the 



gap amplitude as seen in the g2SC case would not be of 
interest here, and partly because we intend to keep an 
interesting resemblance of the present system with elec- 
tric charge turned off and fii introduced to atomic Fermi 
gases with population imbalance created in low temper- 
ature experiments. Note that even for charged systems, 
the gapless superfluid state is generally unstable with re- 
spect to inhomogeneous fluctuations in the order param- 
eter fl7j . 

The typical behavior of the diquark condensate as a 
function of fii is plotted in Fig. [2] for the weak and inter- 
mediate coupling cases. Both are the results at /i q = 0.5A 
and thus the superfluid solution stays at A = 0.05A for 
Parameter I and at A = 0.5A for Parameter II in the 
figure. There appear two branches in the solutions to 
Eq. (fT3"|) for a finite range of /zi, namely, the gapless so- 
lution and the superfluid solution. When these two so- 
lutions meet at a certain /xt, the gap equations cease to 
have any solution. It happens at /ii = 0.099A for Param- 
eter I and fii = 0.99 A for Parameter II in Fig. [2l As we 
shall see from the energy landscape, the thermodynamic 
potential f2 has an inflection point when one solution 
corresponding to a local minimum of the potential meets 
the other corresponding to a local maximum. This means 
that there must exist another state which takes over the 
ground state before fii reaches the meeting point. In fact, 
a first-order phase transition to normal matter is found 
at fij = 0.0355A or, equivalently, 

fii = 0.711A (for Parameter I) , (16) 

where Ao = 0.05A at weak coupling, and [i\ = 0.344A 
or, equivalently, 

\i\ = 0.688A (for Parameter II) , (17) 
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where An = 0.5A at intermediate coupling. As discussed 
in Ref. [24j . this Hi of the first-order transition approxi- 
mately gives the lower bound 5fi2 of the LOFF favored 
region if the upper bound J/n of the LOFF solution, as 
will be discussed later, is greater than 5fi2- The point 
is that the energy gain in the LOFF phase relative to 
normal quark matter is tiny as compared with the scale 
of the change in the energy difference between the su- 
perfluid and normal phases with increasing fii. In order 
for the LOFF phase to be most favorable, therefore, it is 
necessary that the energy gain in the superfluid phase be 
very close to zero, which is only possible when /ii is very 
close to the first-order phase transition point. We note 
that our weak coupling value of Sfi2 = 0.711Ao precisely 
agrees with the known result in Ref. ■ 



C. LOFF Phase 

Let us now move on to the central part of this paper 
that addresses the LOFF favored region in // q -/ii space. 
Here we shall discuss the weak and intermediate coupling 
cases separately. 

In the presence of nonzero q, the free energy has an 
unphysical term proportional to —A 2 q 2 which must be 
subtracted. It is a nontrivial problem how to renormalizc 
this spurious contribution properly in a field-theoretical 
procedure. [See discussions in Refs 0, [l(| about this 
problem in evaluating the Meissner mass.] Our prefer- 
ence is to follow a practical prescription here. After we 
solve the gap equations with a given q to get a = <r(A, q), 
which is not affected by the term ~ — A 2 q 2 in question, 
we can define the subtracted free energy as 

O s (A, q) = r>[a(A, g), A(g), g] - 0[a(A, g), 0, g] , (18) 

in order to determine the gap A and the optimal value 
of q from its global minimum. This prescription is the 
simplest choice consistent with the fact that Sl s should be 
flat in the g-direction in the absence of finite A. The sec- 
ond term in the right-hand side could be — J7[ct(0, g), 0, g], 
but it would make only a negligible difference in numer- 
ical outputs. 



1. Weak coupling case 

We can solve the gap equation <9£1 S /<9A = to obtain 
A = A(g) as a function of q. Figure [3] shows the numeri- 
cal results for A(g) with Parameter I (weak coupling) for 
several values of [i\ in the vicinity of the critical value 
jtii = Sfi2- We can observe from the curve labeled with 
fii = 0.033A that A(g) is smoothly connected up to the 
value at q = which corresponds to the diquark conden- 
sate in the superfluid phase. It is clear from fl s plotted 
in Fig. S| that, when fii = 0.033A, the superfluid phase 
at q = has a much smaller energy than the metastable 
LOFF solution at q ^ 0. Interestingly the LOFF solution 




Pair Momentum q [A] 



FIG. 3: Diquark condensate at weak coupling as a function 
of the pair momentum q, which is calculated for fix — 0.033 A 
(top), 0.034A (middle), and 0.0355A (bottom). 




1 °0 0.05 

Pair Momentum q [A] 

FIG. 4: Energy difference (|18[l at weak coupling as a function 
of the pair momentum q. The three curves correspond to 
those displayed in Fig. [3] 

becomes disconnected from the superfluid solution as /ii 
grows up. Nevertheless, the superfluid phase continues 
to exist separately and is more favorable until the LOFF 
phase turns to the ground state at /ii > 5fi2- 

Then, there occurs a crucial question: How large is 
the upper bound Sfii of /ii at which the LOFF solution 
disappears? The g-range of nonzero A(g) shrinks as fii 
gets larger as indicated in Fig. [3l We find that it is 
located at m — 0.0377A or, equivalently, 

Sfn = 0.755A , (19) 
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FIG. 5: The LOFF favored region bounded by Sfii and 8fi2, 
calculated at weak coupling as a function of /i q . The open 
circles represent the calculated values, whereas the curves are 
interpolations between the neighboring circles by splines. 



where Ao = 0.05A, and that the optimal pair momentum 
at hi = Shi reads q = 0.045A or, equivalently, 



q ~ 1.2(5/ii 



(20) 



Our is slightly bigger than the value reported in 



Ref. [24j. We expect that the discrepancy could be re- 



duced if we choose smaller Ao and smaller Hq- 

It would be intriguing to turn our focus toward a even 
higher density region where the saturation effect is rele- 
vant. Then, the above result for the LOFF favored re- 
gion, 6 Hi — 0.755A > hi > Sh2 = 0.711A , may well be 
altered, and indeed the LOFF window enlarges a bit. 

Figure [5] demonstrates that the LOFF window varies 
with changing quark chemical potential Hq- The shape 
of the curves reflects that of Ao(/i q ) given in Fig.[TJ The 
LOFF window becomes slightly wider in the region where 
the saturation effect is sufficient to make the curves have 
a negative slope. For instance, at Hq = 0.95A the LOFF 
favored region is bounded by 



Shi = 0.080A = 0.73A (> q = 0.95A) , 
Sh2 = 0.069A = 0.62A (> q = 0.95A) , 



(21) 



where Ao = 0.1 10 A. The interval between them is not 
substantial in the weak coupling case, while the behavior 
drastically changes at intermediate coupling as we will 
see shortly. 



2. Intermediate coupling case 

In the intermediate coupling case with Parameter II, 
the gross behavior of the diquark condensate and the 



FIG. 6: Same as Fig.[5]except that calculations are performed 
at intermediate coupling. 



subtracted energy as a function of q is just similar to 
what we have shown in Figs. [3] and |U We thus find 
that such behavior has only a weak dependence on the 
parameter choice. For the LOFF favored region, however, 
the parameter dependence is important as can be seen 
from comparison between Figs. [5]and[5J 

The first-order transition from the superfluid to the 
normal phase occurs above the upper bound of the LOFF 
solution for Hq S 0.7A. This means that the LOFF win- 
dow shuts there and the superfluid phase remains the 
ground state even for hi > SfJ-i until normal quark mat- 
ter appears. This behavior, which is a contrast to the case 
with Parameter I, is caused by stronger coupling rather 
than by heavier mo. In fact, the chiral condensate is 
much smaller than the diquark condensate at /i q ~ 0.7A 
as shown in Fig. [1] 

The LOFF window opens at Hq ^ 0.7A where we find 
the same tendency as the weak coupling case that the 
LOFF window widens by the saturation effect. At Hq = 
0.8A, for instance, the bounds of the LOFF window are 
given by 



Shi = 0.33A ; 
Sh2 = 0.29A : 



0.64A (/v 
0.57A (M q = 



= 0.8A) 

:0.8A) 



(22) 



where A = 0.51A at Hq = 0.8 A. 

We note that our plane-wave LOFF ansatz is only 
the simplest one. If we consider more general three- 
dimensional crystalline structure, the LOFF favored re- 
gion has to be even wider. We would stress that our 
results shown in Fig. [5] strongly support the possibility 
of probing the generic crystalline color superconducting 
phase, along with the isotropic superfluid phase, with 
numerical approach on lattice. In fact, the values of hi 
at which the LOFF phase occurs are less than a half of 
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<7[A]°- 3 .4 



FIG. 7: Subtracted free energy as a function of A and q with 
M = M(A, q) which is the solution to the gap equation. A 
darker colored region has a smaller energy. The LOFF phase 
is the global minimum of the energy. 



(iq, while being larger than the inverse of a typical lat- 
tice size. Possible simulation to detect the LOFF phase 
would not require weaker coupling, finer lattice, or larger 
lattice. 



D. Energy Landscape 

We will take a further look at the LOFF phase found 
at intermediate coupling. Let us pick up a point of 
/i q = 0.8A and fii — 0.3A inside the LOFF favored re- 
gion illustrated in Fig. [6] In order to elucidate how the 
unstable gapless phase is connected to the LOFF phase, 
we picture the energy landscape, namely, we plot f2 s as 
a function of A and q as shown in Fig. [7] 

It is obvious from Fig. [7J that the LOFF phase sits 
certainly at the global minimum of the potential energy 
landscape. The local minimum at A = q = is normal 
quark matter, another local minimum at A ^ and q = 
is the metastable superfluid phase, and the in-between 
local maximum along the A-axis corresponds to the un- 
stable gapless phase. The gapless phase is unstable with 
respect to fluctuations in any direction in the space of q 
and A. Specifically, the instability along the A direction 
is the Sarma instability originally noticed in Ref. [J]. In 
the g2SC or gCFL phase, the remedy against the Sarma 
instability comes from the electric and color neutrality 
conditions. Even with neutrality imposed, however, the 
instability in the q direction is still a problem, which is 
nothing but the chromomagnetic instability in the g2SC 
or gCFL phase. For a neutral superfluid of interest here, 
it is the phase instability. 

Figure [5] is the contour density plot of the numerical 
data in Fig. [7J which allows us to describe graphically 



0.1 



0.2 



0.3 

q [A] 



0.4 



0.5 



0.6 



FIG. 8: Contour density plot of Q. s , which indicate two min- 
ima corresponding to the metastable superfluid phase and the 
ground state of the LOFF phase and one maximum corre- 
sponding to the unstable gapless phase. 



what we mentioned above. The gapless phase has in- 
stabilities leading to normal quark matter, the isotropic 
superfluid phase, and the LOFF phase, among which the 
LOFF phase has the lowest energy. This contour plot 
presumably captures the essence of the instability prob- 
lem and its implications for QCD quark matter. We 
would anticipate that the energy landscape around the 
g2SC phase is more or less similar to Fig. [8] except that 
the gapless phase is stabilized against fluctuations in A, 
that is, the gapless state is a saddle point on the energy 
landscape. It is hard, however, to imagine what would 
happen in the case of the gCFL phase with three fla- 
vors, because there are three predominant diquark con- 
densates. In large dimensional space spanned by more 
variational parameters, it is not straightforward to see 
whether the instability in the gapless phase is directly 
connected to the LOFF phase. There could be several 
distinct LOFF solutions as conjectured in Ref. [l6j]. It 
would be a challenging problem to draw the energy land- 
scape in the three-flavor case and to confirm that the 
instability is responsible for the LOFF phase, as we have 
done successfully in two-color and two-flavor QCD in this 
work. 



IV. CONCLUSIONS 

We investigated the LOFF favored window of fi\ for 
two-color and two-flavor QCD using a mean-field model. 
First, we satisfactorily reproduced the known results at 
weak coupling. We then performed our calculation for 
a set of the model parameters in which the coupling 
strength is intermediate and thus relevant to currently 
available lattice simulations. We found out a nonzero 
interval of fi\ where the single plane-wave LOFF phase 
is energetically more favorable than the isotropic super- 
fluid phase and normal quark matter. For intermediate 
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coupling, we took a close look at the energy landscape 
to see the relation of the Sarma and phase instabilities 
associated with the gapless solution of the gap equations 
with the LOFF phase. 

In this work we do not take account of any possibility 
of the mixed phase or the phas e sep aration for a nonzero 
range of as discussed in Ref. [17| because we turn elec- 
tric charge off in the present system. The chargeless limit 
would be advantageous to lattice two-color QCD simu- 
lation. In the absence of the fermion sign problem, fur- 
thermore, we believe that lattice approach to the LOFF 
phase is feasible [iij . 
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